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r-{ • We investigate a possible asymmetry in the statistical properties of the cosmic microwave back- 

ground temperature field and to do so we construct an estimator aiming at detecting a dipolar 
modulation. Such a modulation is found to induce correlations between multipoles with A£ = 1. 
Applying this estimator, to the V and W bands of the WMAP data, we found a significant detection 
in the V band. We argue however that foregrounds and in particular point sources are the origin of 
this signal. 
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'nI" ■ The Wilkinson Microwave Anisotropy Probe (WMAP) data [l|, Q have raised a number of interrogations concern- 

ing the statistical properties of the temperature field. While these data globally confirm the standard inflationary 

r^ , paradigm '3f| and the concordance cosmological model, they exhibit some intriguing anomalies, particularly concerning 

Qlh the large angular scales. In particular, a huge activity has been devoted to the study of the low value of the quadrupole 
and octopole jj, |M S Q as well as their alignment [a, 13 > two effects that appear to be inconsistent with the standard 

^ • cosmological model. 

j/2 I Besides, many authors have tried to test the statistical properties of the temperature field using various methods. 

d • For instance, it was investigated whether the coefficients agm of the development of the temperature field on spherical 

^ ' harmonics were independent and Gaussian distributed. While, as expected from the standard inflationary picture, a 
' "^ I x^ deviation from Gaussianity seems to be well constrain ed |l C|. there have been some claims that the distribution 
^\ ' may not be isotropic Q llll . Ha . \13L [Ij, [l^ or Gaussian |l6l llTl llsL Il9l | . No real physical understanding of these 

3 ■ measurements have been proposed yet and the origin of these pos sible features is still unknown. Some authors have 



measurements nave been proposed yet ana tne origin or tnese possible features is still unKnown. some autnors nave 
argued in favor of systematic effects [13 while it was argue [13, ^^ that foreground contamination may play an 
important role in these conclusions. 

From a theoretical point of view, there are many reasons to look for (and/or constrain) a departure from Gaussianity 
and/or isotropy of the CMB temperature field. Mode correlation can be linked to non-Gaussianity, in particular due 
to finite size effects |2l|, 122, 123 or to the existence of some non-trivial topology of the universe |25| • While in the 
latter case, one expects to have a complex correlation matrix of the aim, the former leads generically to a dipolar 
modulation of the CMB field ^4]. Such a modulation induces in particular correlations between adjacent multipoles 
(A^ — 1). Similar correlations but with Ai — 2 may also be induced by a primordial magnetic field .26j . In each 
case, the physical model and its predictions indicate the type of correlations to look for and will drive the design of 
an adapted estimator. 
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Inves tigati on of the correlation properties of aim is thus important to correctly interpret previous observational 
results 0, EB 113. Edi EIj Il3 ■ Two approaches are thus possible. Either one defines some general estimators and study 
whether they agree with a Gaussian and isotropic distribution (top-down approach) or one sticks to a class of physical 
models and construct an adapted estimator (bottom-up approach). In this article, we follow the second route and 
focus to the task of constraining a possible dipolar modulations of the CMB temperature field, that is correlations 
between multipoles with A£ = 1. 

In Section ^ we start by some general considerations on the form of the correlation arising from a dipolar mod- 
ulation. We then built an estimator, in Section IlIII adapted to these types of correlations. In particular, we cannot 
use full-sky data and we will need to cut out some part of the sky. The effect of such a mask on the correlations will 
have to be taken into account and included in the construction of the estimator. We apply this estimator to the V 
and W bands of the WMAP data in Section IIVI The V band exhibits an apparent detection. The interpretation of 
this result will require us to compare various masks, and in particular to investigate the effect of point sources on the 
signal to conclude that they are most likely its cause. 

II. GENERAL CONSIDERATIONS 

As explained in the introduction, we focus on a possible dipolar modulation of the CMB signal. Thus, we assume 
that the observed temperature field can be modelled as 

e°'^^(7) = e(7) [1 + e-iyi^_i(7) + eon,o(7) + £+iYi.+im (1) 

where 6 is the genuine statistically isotropic field and where {s-i,eo,£+i) are three unknown parameters that char- 
acterizes the direction of the modulation. The modulation has to be real so that so is real and e+i = — e-i* = e. 
As usual, we decompose the temperature fluctuation in spherical harmonics as 

oo I 

Q(T^) = E E ''(rnYtrnil)- (2) 

1=2 m=-l 

The coefficients aim are thus given by 

aim^ j d^lQ{l)Yim{l)- (3) 

6°'^'' and a°^ are defined and related in the same way. Since 6 is supposed to be the primordial, Gaussian and 
statistically isotropic, temperature field, its correlation matrix reduces to 

{aima*i,m') = CiSii'Smm'- (4) 

A modulation of the form (^ implies that the coefhcients a°^ develop correlations between multipoles with A£ = 1. 
Let us illustrate the origin of this correlation. From Eqs. iQ and Q, we deduce that 
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^ J2 «^'-' E ^W <^^in,m{l)Y',n'{7)Yi4j). (5) 

i'=2m' = -i' i=-l ■^ 

The integral can be easily computed by using the Gaunt formula [see Eq. (|A7ll ] to get 

«t^ = aim + a/JE ^^("1)" E «^M V{2i+l)i2L + l) ( _i M 1) ( J ) • (6) 

Because of the triangular inequality, the Wigner 3j'-symbols are non zero only when L ^ £ ± 1 and AI — m ~ i so 
that all^ is in fact a sum involving aim and ai±im-i- It follows that it will develop £—{£+!) correlations that can 
be characterized by the two quantities 

^{m, = {aim ai+im) ' (''') 

-^£m = {aim ai+im+l) (8) 



which will be non zero respectively as soon as Eq or e are non zero. Using the expression © and the property Q of 
the primordial field, we deduce that 
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Interestingly, these forms indicate how to sum the Dgm in order to construct an estimator. This construction will be 
detailed in the following section. 

III. MATHEMATICAL CONSTRUCTION OF THE ESTIMATOR 

The previous analysis is illustrative but not suitable to be applied on real data. In particular these data will not 
be full sky and we have to take into account the effect of a mask (see e.g. Ref. |23). Such a mask, that arises in 
particular because of the galactic cut, will induce correlations in the coefficients afj^ that are described in § IIII Al 
We design the mask in order to protect the correlations that originate from the modulation fS IIII Bl and !^ IIII C|) and 
finish by presenting the construction of our estimator in the most general case (§ IIIID|I . 

A. Mask effects 

The temperature field is observed only on a fraction of the sky. We thus have to mask part of the map so that the 
temperature field is in fact given by 



e°^^(7) = e(7) 



1+ ^£.ri..(7) 



W{j) (11) 



where 14^(7) is a window function, referred to as mask, indicating which part of the sky has been cut. We decompose 
1^(7) in spherical harmonics as 

M^(7)=I]w^^ml£m(7)- (12) 

em 
W{j) being a real valued function, it implies that w^^^ = (— l)™u;£_m. We deduce from Eqs. Hll|) and ||2Jl that 

a,t^=S,„ + ^e,4:l (13) 

i 

where agm are the coefficients of the masked primordial temperature field 0(7) — 6(7)1^(7), 

Y^ a^imi Yl '^^2m2 / d^7y«imi(7)yf2m2(7)y/m(7) (14) 

iimi ^2"i2 

and the effects of the modulation are encoded in the correction 

^£m = Yl "'^i"! Yl ^^2^2 / <i^iYi,md7)ye2m2{7)yii{l)Yi*,ri{l)- (15) 



aim — 

timx ^2"l2 



(?2"l2 



Interestingly, aim can be shown to be obtained from aim by the action of a Kernel Kf, 



= Y <^i^mXr'- (16) 
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This kernel is defined by 

and can be explicitly computed by using the integral ljA7|l to obtain 



^^m - (-1) 2. ^^.™.V ^^ I mi m2 -m ) I ) ' ^^^^ 

The contribution arising from the modulation can be computed by using the integral l)A7|l to get 



4:i - ^i^i-^r E «-^ V21TTV2LTT [M--m[ooo- (i«) 

One can check that the relation Q obtained without taking into account the effects of the mask still holds if one 
replaces aim by a^m- The complications arise from the fact that Hirn does not satisfy the property l@J because of the 
action p(j|) of the Kernel. 

B. Choice of the mask and properties of the masked quantities 

We now need to specify the form of the mask. First, let us note that when W{^) =constant for all 7 then one 
trivially recovers that Hgm = o-tm because W — uioo^oo so that 
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Since we are looking for £ — (€ + 1) correlations, we would like to design a mask that does not involve the same 
correlations for 'dim and that is not rn-dependent. A solution is to impose that ^^(7) is a function of 9 only and that 
it is north-south symmetric, that is 

w{^) = w{e), w{tt -0)^ w(e). (19) 

Since Y^ol^r — 9) = (— l)^Y{o(^)i these conditions imply that 

W{i) = E^^^^^^T^^'^"^^^ ^ E^^^^^^o(^) (20) 

The simplest example of such a mask is obtained by considering a function which is constant and vanishes on an 
equatorial strip of latitude 9c G [0,7r/2]. This implies that the multipolcs of the mask are given by 



wq = v47r/Xc, (21) 

m = y^^^i^-^[P.-i(/ic)-P.+i(A^c)] (22) 

where fj,c = cos9c- In particular, it can be seen that when 9c -^ 0, that is when the size of the mask vanishes, this 
mask satisfies wg — > \/4n6io when /ic — > 1. The function wi is depicted on figure ^ for galactic cuts of 10, 20 and 30 
degrees. 

The results derived in the following sections are not dependent on the particular choice of the mask as long as it 
satisfies the symmetries (I19|l which ensure that the coefficients of the mask do not depend on m and vanish for £ odd 
(see Eq. 1201 ■ 

C. Properties of the aim 

Whatever the choice of the mask, as long as it satisfies the properties (|19f) , the general expression of the coefficients 
aern of the decomposition of Q, are given by 



de,n = /.kya.™ + (-1)'"V2ZTTE V ^%^"^^'» E V^^^+^^i. ( t ^0 -m ) ( ^0 t) ) ^^S) 

il ^ l2^0 ^ ^ ^ ^ 
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FIG. 1: (left) The coefficients of the decomposition of the mask function on spherical harmonics for different galactic cuts, 
/ic ~ 10, 20 and 30 deg. (respectively black, red and blue lines). Only the non- vanishing, that is even, multipoles are plotted, 
(right) The ratio we/wo decreases rapidly and becomes typically smaller than 0.1% for £ larger than 20 (plain=10 deg., dot=20 
deg., dash-dot=30 deg.). 



with /sky = Wo/ v47r is the fraction of the sky that is covered. From this expression, we deduce that their 2-point 
function is given by 



{aima},„,,) = (5,w <^ Q/iyfe' + /sky \Q{ll, t', Tfi)Ct + e(/, ^, ni)Ci,\ + ^ C,,g(4, ^, m)g{i^,i' , m) 



where the function Q{£, £' , m) is defined by 



gie,£',m) = (-i)'V(2^ + i)(2^' + i)E 
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(24) 



(25) 



It follows from Eq. (|24|l that there is no 771-coupling arising from the mask (because it has no azimuthal dependence) 
and we can define the correlation matrix of the masked temperature field as 



The angular power spectrum of the mask field is then defined as 
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and is explicitly given in terms of the primordial angular power spectrum by 
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Let us now turn to the i-{i + 1) correlators. The first term in Eq. (|24|l vanishes. Then, one can check that 
Q{i^ £ + l,m) vanishes because the triangular relation of the Wigner-3j symbols implies that £2 = ±1 but for odd £2 
wi^ vanish. To finish, the contribution of Q{£i^i,m)Q{£i,i + 1,771) in the sum also vanishes because £2 is even and 
the sums ^1+^2+^+1 and £1 -I- ^2 + ^ have to be both even, which is impossible. In conclusion 



{ahna*i+im') - 0- 



(29) 



As expected from our construction, the mask does not generate £-{£ + 1) correlations. 

To finish, let us stress that the mask will induce some £-£ + 2 correlations that can be characterized by introducing 
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Indeed, when W = Id, 2C1 — 0. 



D. General construction 

Starting from the relation ifT^ and the expression l(TS|l . we deduce that the two quantities defined in Eqs. H7I8|) 
generalize to 
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(31) 



with i = 0, 1, when the mask effects are taken into account. This expression is defined for m = —£ . . .£ even if 
2Ci-im+i is not defined for m — £ and m = £ — 1 and 2Ci-im for m = £ because the Wigner-3j symbols that 
multiply these terms strictly vanish. From this expression, we define 
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(32) 



Now, it can be checked, after some algebra, that 
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where the quantities C/ , ^Cg and 2C1 have been defined by 



c« = 4. E v/^TT) 



iV, 



(i) 



^ + 1 1 £ 
m i —m — i 



2^ 



{i)A _ 



2C, 



{^)B 






1 ^ 



Clm + C^+lm+i 



'+2 1 ^+1 

TO i —m — i 



iCtm 



£ I £+1 



£-\ \ £ 



m I —TO. — I I \ m + t —I —TO 



2Ctm 



(34) 
(35) 
(36) 



with the coefficients iV^ given by 
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It follows from these results that we can consider the estimator 



£;« = 



1 ^ 



£ 1 £+1 
TO i — ?7i — i 



(-1) 



i.+7n+i obs obs* 



(38) 



that satisfies by construction 



I^)^Df^ 



(39) 



We will apply this estimator to the WMAP data in the following sections. 
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FIG. 2: Df measured on the WMAP data (W band) [green triangles]. Blue triangles are the measurements on a simulated map 
with £o = 0.2 and the red hnes are la error bars. 



IV. DATA ANALYSIS 



The proposed estimators have been implemented numerically, using the Healpix ^52| package for the pixelization 
and the fast spherical harmonics transforms, and applied to the co-added data of the WMAP V and W bands (resp. 
70 GHz and 90 GHz) where most of the signal is of cosmological origin. We implemented the estimators as described 
by equations l|5^ to (|S5|l . 

The quantities Cj ' \ ? Ct and 2C/ have been computed using the best fit LCDM theoretical power spectrum 
of the WMAP data [23, and were not computed on the data itself to avoid ratios of random variables. To assess the 
statistical significance of the measured values of Si, we made 1000 simulations of WMAP data in each of the V and 
W bands according to a sky model with no modulation. 

The results of the analysis of the V and W bands are summarized on Figs. [21 to [S] Fig. [3 depicts the measurement 
of Dg on the W band. We sum this measurement of two bands of £ (respectively 20-100 and 100-300) and compare 
with 1000 simulated WMAP data. We perform the same tests on the V band (Fig.2]and[Sl). The apparent detection 
in the V band without clear counterpart in the W band suggest a non cosmological contamination. Determining its 
origin requires to performed more tests. 

This contamination can be a priori from two possible sources, Galactic or extragalactic. To check if the correlations 
detected in the V band are of Galactic origin, we apply the same estimator to the half sum and half difference of the 
V and W bands, that is 



^_ w + v w-v 

o — ~ , J-J — ~ 



(40) 



and repeat the whole procedure on 1000 simulation in each case, where the simulations contain only CMB and noise 
according to the WMAP specifications. The advantage of the half difference of the bands is that it should (up to 
calibration errors) eliminate the CMB signal completely at large scales, hence eliminate the main source of variance 
at these same scales, where the Galactic signals are expected to dominate. Indeed, the power spectra of Galactic 
emissions usually scale as Ce oc i~°', with 2 < a < 3 (see e.g. Ref. z^). The half sum results, summarized in table I, 
are in between those of the V and W bands, which is coherent with the assumption of the detection being caused by 
a foreground source of electromagnetic spectrum different from the CMB fluctuations. 

More importantly, the half difference results do not show a strong correlation detection at large angular scales, 
in contradiction with the assumption of the Galactic foreground contamination being the source of the detected 
correlations in the V band. 

However, this half difference test does not work that well if the source contaminants are of extragalactic origin, 
since the power spectra of extragalactic foregrounds resemble that of the noise. In this case, the contamination is 
expected to increase with increasing multipole number, which seems to be the case for the V band (see Figs. 0] and 

The difficulty of extragalactic point sources contamination is that these sources (quasars and active radio-galaxies) 
are distributed more or less uniformly across the sky, which renders their masking by an azimuthally symmetric sky 
cut impossible. However, the WMAP team provides with their data sets "taylor cuts" that blank out the resolved 
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FIG. 3: Comparison of 1000 simulations with the WMAP data. We use the W band and sum the multipole between ^ = 20 
and £ = 100 (top panel) and between i = 100 and i = 300 (bottom panel). The color code is identical to the one of Fig.|5|that 
is the green line correspond to the measurement on the WMAP data and the blue line the measurement on a simulation. 



point sources of largest flux. Of course, the dipolar modulation estimators designed in the preceding sections do 
not apply stricto sensu to these arbitrary masks, but one can hope, given the small fraction of sky removed at high 
latitude, that the broken symmetry of the mask will be a small perturbation in the computation of the e's, so that 
the estimators keep their general validity, up to a possible small bias (see Fig. |S1 for a comparison of the coefficients 
We of the two masks). 

This assumption can be checked on a simulated sky with a known dipolar modulation, where a WMAP V-band 
noise is added to the signal, together with the taylor mask. We chose the most conservative mask provided by the 
WMAP team (kpO), and applied it to a simulated sky of known modulation {eq = 0.2) as described above, together 
with the V-band data. The results are shown in figures IHl and [7\ Again, the estimators have been applied to 1000 
simulations of the V band with no dipolar modulation, with the same kpO mask applied, to estimate the statistics of 
the V-band data results. 

Several observations can be made at this point. First, comparing these results with those obtained in the V-band, 
but with the azimuthally symmetric 20° cut (figure 0J), one can check that the estimators give very compatible results 
for the simulated dipolar modulation (blue squares) . This comforts our assumption that changing the cut sky to the 
kpO mask is a small perturbation for the modulation estimators. 

Secondly, comparing the same figures but this time looking at the data (green squares), one can see that in the case 
of the 20° cut there is a large trend at high fs in ei that disappears when the kpO cut is used. This is confirmed by 
the results of table U where one can check that the tentative detections of in the V band using the simple cut become 
statistically insignificant when using the kpG cut. 



V. DISCUSSION AND CONCLUSIONS 

In this article we have proposed an estimator designed to detect a possible modulation of the CMB temperature 
field, or equivalently A^ — 1 correlations. The effects of cutting part of the sky were discussed in details and we 
applied this estimator to the V and W bands of the WMAP data. 

The results of our analysis are summarized in table ^ which gives the amplitude of the modulation coefficients on 
the WMAP data and a corresponding test case with eg = 0.2, Re(£i) = Im(ei) = 0. All values are given in standard 
deviation units, estimated on 1000 (signal+noise) simulations in each case, with no modulations. 

While the V band seems to exhibit a marginal detection, further tests such as the study of the half sum and 
difference of the two bands and the effect of point sources have led us to conclude that this detection should be 
inferred to the effect of point sources contamination. In this analysis we have used the kpO mask which does not 
satisfied the symmetries of the mask required for our estimator to be unbiased. Nevertheless, our estimator seems to 
be well suited for the analysis, even with the kpO mask. 

To back up this interpretation we have performed two last tests. First we added to a simulated CMB map without 
modulation and with noise the 208 sources resolved by the WMAP experiment and then smoothed with the correct 
beam. Second we added to the same simulation the 700 circular region that are cut in the analysis of the V band 
in the WMAP analysis. Both simulations, while analyzed as the previous data with an azimuthal mask of 20 deg., 
exhibit an excess of signal for Eq and Im(ei) in the same range of multipoles than obtained on the analysis of the V 
and W bands (Figs. El and 0}. Interestingly, the signal of Re(£i) is not affected and is identical to the one of Figs. El 
and El Indeed, the signals have not exactly the same amplitude as the ones obtained from the analysis of the V band 
but they exhibit the same trend on the the same scales. Also, it has to be stressed that with a cut of 20 deg. the Large 
Magellannic Cloud (galactic latitude of 20 deg. and more and longitude of deg.) and a part of the 112 Ophucius 
region should contribute and that we have not included them in the simulations. This could have enhance the signal. 

In conclusion, the set of analysis performed in our study tend to show that the A£ = 1 correlations that appeared 
in the analysis of the V and W bands of the WMAP data are due to foreground contaminations and most likely by 
point sources. The direction of the detected modulation will, in that interpretation, characterize the anisotropy of 
the distribution of these sources. 
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APPENDIX A: INTEGRALS OVER SPHERICAL HARMONICS 



We have evaluated integrals over n spherical harmonics (see Ref. 31j). When n = 1 or 2, these integrals are trivial 

d^7Y«m = VinSioSmO (Al) 
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FIG. 4: Di measured on the WMAP data (V band) [green traingles]. Blue triangles are the measurements on a simulated map 
with Eo = 0.2 and the red lines are Icr error bars. 
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To go further, one solution is to use the decomposition of the product of 2 spherical harmonics as 

>£imi(7)^<?2m2(7) = 2J V 4 /2L + 1) '^ii0£20C|;mi£2m2^iM(7) 
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where the Cf ^ » are the Clebsch-Gordan coefficients that can be expressed in terms of Wigner 3j symbols as 
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It is easy to generalize Eq. (IA3(I to a product of n spherical harmonics 
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We deduce, using Eq. IjASp and the integral (|A2|) that 
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FIG. 5: Comparison of 1000 simulations with the WMAP data. We use the V band and sum the multipole between ^ = 20 
and i = 100 (top panel) and between i = 100 and i = 300 (bottom panel). 
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FIG. 6: Di measured on the WMAP data (V band) [green triangles], using the taylor mask kpO to blank the main point 
sources. Blue triangles are the measurements on a simulated map with eo = 0.2 and the red lines are la error bars. 
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TABLE I: Summary of the data analysis performed in this article. It concerns the two bands V and W, their half sum (S) and 
half difference (D). The band V has been analyzed with two masks to emphasize the effect of the point sources on the result. 
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FIG. 8: The coefficients of of the decomposition of the mask function on spherical harmonics for different galactic fic ~ 10, 20 
and 30 deg. (respectively black, red and blue lines) compared with the ones of the kpO mask (green) . 
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FIG. 9: The analysis of the simulated maps in which (top) the 208 resolved sources of the WMAP catalog have been added 
and (bottom) where the 700 point sources (resolved and unresolved) of the V band have been added. In both cases, the signal 
of Re(ei is not affected while the signal of £o and Im(ei) exhibit patterns that are similar to the ones obtained in our analysis 
of the V and W band on the same angular scales. 



